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The revival structure of wave packets is examined for quantum systems 
having energies that depend on two nondegenerate quantum numbers. 
For such systems, the evolution of the wave packet is controlled by two 
classical periods and three revival times. These wave packets exhibit 
quantum beats in the initial motion as well as new types of long-term 
revivals. The issue of whether fractional revivals can form is addressed. 
We present an analytical proof showing that at certain times equal to 
rational fractions of the revival times the wave packet can reform as a 
sum of subsidiary waves and that both conventional and new types of 
fractional revivals can occur. 
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Among the distinctive quantum features of wave-packet evolution in a potential is 
the formation of fractional revivals. These appear when the packet reforms as a sum 
of distinct subsidiary waves some time after its initial collapse. Fractional revivals 
have been observed experimentally in the motion of atomic Rydberg wave packets 

0, |3| and in molecular vibrational wave packets |4jj. In these systems, the wave 
packet is formed via excitation with a short laser pulse, and its motion is detected 
using a pump-probe configuration Q. The revival structure of a free hydrogenic 
wave packet is determined by the dependence of the component-eigenstate energies 
on the principal quantum number. Similarly, the vibrational states of a molecular 
wave packet depend only on the vibrational quantum number. 

For the case of quantum systems with energies depending only on one quantum 
number n, a generic proof has been given showing that fractional revivals form ||. 
The proof applies to packets that are a superposition of states strongly peaked around 
a central value n of n. Initially, the motion of a wave packet of this type is periodic 
with period T c \ equal to the period of the corresponding classical motion. After several 
classical periods, the wave packet spreads and collapses. The subsequent formation 
of revivals is characterized by a time scale t Tev called the revival time. Fractional 
revivals occur at times that are rational fractions of t rev . The associated subsidiary 
wave packets exhibit a periodicity in their motion equal to a fraction of T c \. 

A full revival can be viewed as a special fractional revival occurring when the wave 
packet reforms as a single entity evolving with the classical periodicity. For quantum 
systems with energies that are quadratic in n, the full revivals are perfect in the sense 
that the wave function at time t = t rev is exactly equal to the wave function at t = 
0. For other systems the revivals at t — t rev are imperfect, and new structure emerges 
after several cycles of collapse and revival. A detailed analysis involves introducing a 
superrevival time scale t sr and showing that the interplay between t sr and t rev results 
in a new sequence of fractional and full superrevivals for t ^> t rcv ||. Among these, on 
time scales much larger than the superrevival time scale t ST , are the long-term revivals 
already predicted in early work on Rydberg wave packets that appear when the 
phase in the wave function is an integer multiple of 2n [ |T0|| . 

The above analyses can be applied provided the physics of the wave packet is 
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controlled by a single quantum number. As an example, consider a radial Rydberg 
wave packet in an alkali-metal atom. The energies of the atom have quantum defects 
and hence depend on more than one quantum number. However, for large enough 
principal quantum number and fixed angular-momentum quantum number, the quan- 
tum defects are constant to an excellent approximation. The fractional revival and 
superrevival structure of a wave packet formed as a superposition of, for example, p 
states can therefore be treated using the above methods || taking the quantum-defect 
shifts into account analytically JL1]. A theoretical description of these wave packets 
as squeezed states can be used to extend the treatment to elliptical wave packets in 
alkali-metal atoms |12| . 

To study more complicated cases, such as Rydberg wave packets in external fields, 
an analysis is required that allows for nondegenerate energies depending on more than 
one quantum number. This is the topic of the present paper. 

We focus attention here on the dynamics of a wave packet ^(t) formed as a 
coherent superposition of states 4> nin2 with energies E nin2 depending on two quantum 
numbers n\ and n 2 : 

*(*) = c ni „ 2 0„ in2 exp [-iE nin2 t] . (1) 
711,712 

Here, c nin2 are complex weighting coefficients. Note that n\ and n 2 may represent 
effective quantum numbers in, for instance, a perturbative treatment, for which the 
states 4> nin2 may be only approximate energy eigenstates. Our analysis in the present 
work is therefore relevant for experiments studying the behavior of Rydberg wave 
packets in external fields, where the energy degeneracy in the superposition of states 



is lifted 00 [15], ]I6|, 0. 

We suppose that the superposition is weighted around central values n\ and n 2 
of the two quantum numbers. The weighting coefficients c nin2 are therefore taken to 
be strongly peaked around a central value of the energy E nin2 , as would result from 
an excitation process using short laser pulses. The energy can then be expanded in a 
Taylor series as 

f dE\ , _ N (dE\ 
E nin2 ~ E nin2 + 1^-1 (ni-ni) + I I fa - n 2 ) 
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(ni - ni)(n 2 - n 2 ) + . . . . (2) 

ni,n2 

Since the energy depends on two variables, the expansion contains a mixed derivative 
term at second order. 

When Eq. (Q) is substituted into (|]) and an irrelevant overall phase is disregarded, 
it follows that each derivative term corresponds to a distinct time scale. We define 

(1) _ 2tt (2) _ 2tt 



t (l) = — f(2) = ^ M 

4lv = 7 Q2 E \ > ( 5 ) 

V9ni9n2yni,n2 

where we have assumed for simplicity that the quantum numbers take integer-spaced 
values. For each quantum number, there is a classical period and a revival time scale. 
In addition, the mixed derivative term generates a third time scale t^) , which we call 
the cross-revival time. This term has no analog in the conventional revival structure 
of free Rydberg wave packets. Keeping terms to second order then gives for the 
expression 



*(*) = J2 c klk2 (f) klk2 exp 



. / kit k 2 t k\t k\t kik 2 t^ 

~ 2m { ^ol + Tfm +717 + 727 + 712T 

\ cl cl (rev (rev t-rev / 



(6) 



where k\ = (ni — ni), k 2 = {n 2 — n 2 ), and where the coefficients c klk2 an d eigenstates 
4> kl k 2 are related to the analogous quantities in Eq. flU) by shifts in the subscripts. 

Consider the time evolution of the wave packet. When t is small, the dominant 
terms in the time-dependent phase in Eq. (|6]) are the first two. They produce beating 
between the classical periods T cl and The latter are said to be commensurate 



if 



lf = \lf , (7) 
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where a and b are relatively prime integers. When the classical periods are commen- 
surate, the wave-packet displays a period on short time scales given by 

T d = bT^=aT^ . (8) 

If two classical periods are incommensurate, the initital motion of ^{t) is not exactly 
periodic. 

On larger time scales the three quantities t[H, t^, and t££) play a role in control- 
ling the behavior, eventually inducing the wave-packet collapse. It follows from the 
form of the three second-order terms in the time-dependent phase that full revivals 
appear if the three revival times t^, t^, and are commensurate and satisfy 

f (i) = £ f (2) = £ f (12) (q \ 

"rev d ev y rev ) K^J 

where c, d and e, / are pairs of relatively prime integers. 

When Eq. ([]) is obeyed, there is a time t = t rev at which all three second-order 
terms in the phase become integer multiples of 2ir. Near this time, a full revival 
occurs: the phase is governed once more by the first-order terms, so the wave-packet 
motion and shape resemble those of the initial wave packet. Typically, the full revival 
time t rev is a multiple of tQl, t^l, and t££ '. However, it may equal one or more of the 
three component revival times for particular values of the integers in 

The appearance of fractional revivals requires that the wave function \I/(t) in Eq. 
(|) be expressible as a sum of subsidiary wave functions. This is possible only at times 
t = tfrac that are simultaneously irreducible rational fractions of all three revival time 
scales. We therefore write the fractional-revival time scale as 

fr _ El f (i) _ £a f (2) _ PRM2) (1Q] 

f-frac — ''rev — ^rev — ''rev ■ 

qi <?2 qn 

Here, the pairs of integers (pi,qi), (^2,^2), and (^12,^12) are relatively prime. Note 
that the condition @ for the appearance of a full revival is automatically fulfilled 
if times if rac can be found that satisfy Eq. ([!(]). The formation of full revivals is 
therefore a special case in the analysis below. 

If a set of integers obeying ([T0|) exists, then the issue of the appearance of fractional 
revivals at t = tf rac becomes pertinent. To make further progress, we introduce a 
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doubly-periodic function depending only on the classical time scales: 



Ai(ti,h) = c fcifc 2 0fcifc 2 ex P - 27r ^ — m + 



ki,k 2 I \ J cl J cl 




(11) 



Here, t\ and t 2 are dummy variables. The function ip c \(ti, t 2 ) is periodic in t\ and £2 
with periods and X^ 2 - 1 , respectively: 



The functions V'cit^i,^) have no immediate physical significance. However, at the 
time t = ti = t 2 the function ip c i(t,t) is equal to the wave function ^(t) of Eq. (^) 
when the second-order terms in the phase are omitted. Thus, the functions ip c \(t,t) 
exhibit the beating of and Note that ^(0) = -001(0,0), so for small times 

the behavior of is approximately matched by that of ipci{t,t). 

The functions V'ci^i)^) are useful because we can prove that, near the times tf rac 
given in Eq. (|TUD, the wave packet can be written as a sum of subsidiary waves ip c \ 
with arguments shifted relative to t by certain fractions of the corresponding periods. 
The proof involves examining the cyclic properties in k\ and k 2 of the second-order 
terms in the time-dependent phase of and demonstrating that the functions ip c \ 
with shifted arguments have identical cyclic properties. The result means that the 
ip c \ with shifted arguments form an acceptable basis for the expansion of at the 
times t frac . 

To examine the cyclic properties in k\ and k 2 of the second-order contributions 
to the time-dependent phase in at t = tf rac , we write these contributions as 
exp (-2iri9 klk2 ) , where 



(12) 



(13) 




(14) 



We seek the minimum periods l\ and l 2 such that 



(15) 



Note that the relations 



(i) = HA12) 



(2) = ?X(12) 



(16) 



rev _ rev 
■Si 



rev "rev 
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follow from fllCf ) with 

n _ qiPi2 r2 - q2Pl2 

si Pxqvz ' s 2 p 2 qi 2 
The periods l\ and l 2 must then satisfy 



(17) 



Elll + —hh + ^hk 2 = (modi) , (18) 
qi 9i 

-M f2^2 H t 2 «i = (modi) . (19) 

<?2 <?2 ?2«2 

These relations can be satisfied by choosing l\ = qiS\ and l 2 = 92 ^2- However, smaller 
values might exist if cancellations occur in some ratios. 

The second-order contributions 6k 1 k 2 are cyclic in h\ and k 2 with periods l\ and 
I2, respectively. The functions ip c \(t + siT^/h, t + s 2 T^ // 2 ) with shifted arguments 
have the same periodicities in k\ and k 2 . Here, ip c \ is defined in Eq. ( |TTD and Si = 
0, 1, . . . , l\ — 1 and s 2 = 0, 1, . . . , l 2 — 1. This suggests the set of functions ipci with 
shifted arguments can be used as a basis for an expansion of \I/(t) near the time tf rac . 

Explicitly, the wave packet at the times t pa tf rac can be written as 

f^EE^+r^t+r 1 ?) ■ ( 2 °) 

The integers l\ and Z 2 depend on the values of if rac . They are obtained by solving the 
conditions ( fL8|) and (|19"D. The coefficients a SlS2 are given by 



1 / r, •/> \ /r> - S l \ In - S 2 



11 exp(-27rz6» Klft2 )exp yiixi—K X exp [2^1— k 2 ) , (21) 

( 1'2 Ki =0k 2 =0 \ h / \ l 2 J 

where 9 K1K2 is given by Eq. (0). If the expressions for a SlS2 and the ip c \ are substituted 
into Eq. (^Cj), the expansion reduces to Eq. (|6]) evaluated at the time t = £f rac , which 
completes the proof. 

The expansion in Q2!]| ) shows that the wave packet can be written as a sum of 
subsidiary wave functions at t pa i frac . However, it can be seen from Eq. ( BD[ ) that the 
subsidiary wave functions typically have distinct shifts in the two arguments and hence 
do not follow the initial evolution of ^(t), which is given by ip c \(t,t). Nonetheless, 
if near these times the shifted subsidiary wave functions are macroscopically distinct 
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and localized in space, then ^(t) has a periodicity that is a fraction of the classical 
time periods. 

The detailed behavior of the fractional revivals is governed by the form of the 
expansion coefficients a SlS2 given in Eq. These coefficients depend on the func- 

tion defined in Eq. (|Tip, which in turn depends on the commensurability of 
the three component revival times. An interesting special case occurs if the cross- 
derivative term in the energy expansion vanishes, d 2 E / 'dnidn 2 = 0, whereupon there 
is no cross-revival time tjffi . The coefficients a slS2 can then be written as a product 
a *in a &?i where aW and each separately have the form of expansion coefficients 
for the conventional fractional revivals ||. Nonvanishing values of acquire the 
same norm for all s\. The same holds for af^. This means that for zero cross- revival 
time the coefficients a SlS2 are equal to a product of the corresponding coefficients for 
systems with conventional fractional revivals. Nonetheless, the function ip c \ in Eq. 
(|TT|) cannot be separated into a product of one-dimensional functions unless there is 
no entanglement of the states in the superposition. Only in this special case can the 
fractional revivals for systems with energies that depend on two quantum numbers 
be written as a product of two conventional fractional revivals. We refer to super- 
positions having these special properties as totally separable wave packets. For such 
packets, it follows that if Ni subsidiary wave packets appear for one component while 
N 2 appear for the other, then the fractional revival for the full system consists of 
N1N2 subsidiary wave packets. 

If the cross-revival time is nonzero, the expansion coefficients a SlS2 cannot be 
written as a product a^a^. Their form is more complicated, and typically they have 
different norms for distinct s\ and $2- The fractional revivals are therefore less likely 
to cause distinctive periodicities in the autocorrelation function. However, for certain 
circumstances only a small number of subsidiary wave functions enter the expansion. 
This can happen, for example, if the values of l\ and I2 are small or if some of the 
a SlS2 coefficients vanish. In this case, fractional periodicities should appear in the 
autocorrelation function. 

Typically, the conditions fllOD for the existence of a time tf rac at which fractional 
revivals can occur are quite restrictive. For most values of and t^v , 110 
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solution exists. However, for many systems with energies that depend on two quantum 
numbers and that are of experimental interest, there is at least one tunable parameter 
such as a field strength or a length scale. If the energy depends on the tunable 
parameter, then so in general do the component revival times. It may therefore be 
feasible to tune the system so that a particular fractional revival exists. Since only 
certain commensurabilities are likely to hold, a partial spectrum of fractional revivals 
should be seen. For different values of the tunable parameter, we expect different 
partial spectra of fractional revivals. 

As an example of a totally separable wave packet, consider a superposition of 
energy states for a particle in a two-dimensional box with periodic boundary con- 
ditions. The box is taken to have lengths L\ and L 2 in the x and y directions, 
respectively, and the mass of the particle is chosen to be that of the electron, i.e., 
one in atomic units. The eigenfunctions are separable in x and y, with eigenenergies 
E niri2 = 2n 2 {n\Li 2 + n\L 2 2 ). The wave packet is taken as a superposition of states 
centered on fi\ and n 2 with weighting coefficients c nin2 = c^c^, where |c^| 2 and 
|c^| 2 are both gaussian functions with widths <j\ and a 2 , respectively. 

Since the energy separates, there is no cross-revival time. Furthermore, since more 
than two derivatives of E nin2 vanish identically, there are no superrevivals for this 
system. This means the full revivals are perfect revivals for which the absolute square 
of the autocorrelation function is one 0. The only nonzero revival times are T<jj = 
L\j2'Kfix 1 = L 2 ,/27rn 2 , t^H = Lf/ir, and t^l = L 2 /ti. The commensurability of 
the revival times may be tuned by adjusting the ratio L\jL\. 

For definiteness, consider a superposition with n\ = n 2 = 18 in a system with 
Li = V3/2 and L 2 = 1. The ratio of the revival times is then t^H/t^l = 3/4. We 
expect a full revival at t = t rev = 4t^| = 3t^| = 3/n ~ 0.95 in atomic units. Figure 1 
shows the absolute square of the autocorrelation function |A(t)| 2 = 1(^(0) ^(t))! 2 as 
a function of time. A perfect full revival is evident at t = t TCV . At t = t^H ~ 0.24 a.u. 
and t = 3*W ~ 0.72 a.u., corresponding to the choices pi/qi = 1, = 3/4, and 

P1/Q1 — 3> P2/Q2 = 9/4, respectively, two subsidiary wave functions form. At these 
times \A(t)\ 2 = 1/2, as can be seen in Fig. 1. At t = 2t« ~ 0.48 a.u., corresponding 
to Pi/qi — 2 and p 2 /q 2 = 3/2, one wave function forms in each component. However, 
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this time t = t^l/2, which means the second-component subsidiary wave is a half 
cycle out of phase. Thus, as seen in Fig. 1, the autocorrelation A(t) is approximately 
zero at this time. 

As a second explicit example, set L\ — 1 and L2 = V3, keeping fi\ = n 2 = 18. 
These dimensions produce the ratio 4el/4el = 1/3, and as a result different types 
of fractional revivals will occur. The autocorrelation function for this case is shown 
in Fig. 2. A perfect full revival occurs at t = t vev = 3tQl = 4t^\ = 3/ir ~ 0.95 
in atomic units. Fractional revivals consisting of three subsidiary wave functions 
form at multiples of t^H/2 ~ 0.16 a.u., with the exception of the time 3t^H/2. At 
this time, pi/qi = 3/2 and = 1/2, and each component consists of a single 

wave packet forming one-half cycle out of phase with the classical motion. This does 
not result in a perfect revival because the cycles are not equal, Tjp/2 7^ /2, 
so the two components are out of phase. Additional fractional revivals consisting 
of four subsidiary wave functions are visible at times t = 3^1/4 ~ 0.24 a.u. and 
t = 9tW/4 ~ 0.72 a.u. 

The analysis in this paper is experimentally testable. As an explicit example of a 
physical system for which the energy depends on two quantum numbers and a cross- 
revival time exists, we can consider a Stark wave packet. This type of packet can 
be produced using a static electric field applied to an atom, which shifts and splits 
the energy levels. Applying a short laser pulse then creates a coherent superposition 
of Stark levels. The energies in atomic units for hydrogen and a weak electric field 
are E n k = —l/2n 2 + 3nkF/2, where n is the principal quantum number, k is the 
difference ri\ — 122 between the parabolic quantum numbers, and F is the magnitude 
of the electric-field strength. 

Stark wave packets have unique features because the quantum number k is even 
or odd according to whether n is odd or even and because for fixed n the adjacent k 
values differ by two units. These properties introduce some additional complications 
in the analysis of fractional revivals and produce novel interference patterns in the 



interferograms of Stark wave packets [18|. A second-order Taylor expansion of the 
energy E n k with respect to n and k can still be performed, but some changes occur in 
the definitions of the controlling time scales. Since the quantum number k is either 
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even or odd, a Stark wave packet effectively evolves as two separate wave packets. 
At the fractional revival times, the packet separates into sums over sets of subsidiary 
wave functions consisting of superpositions of even-ra and odd-n states. These wave 
functions evolve with distinct relative phases and obey distinct periodicity relations. 
This results in additional interference patterns in the autocorrelation function, which 
can be tested experimentally. 

One of us (R.B.) would like to thank Colby College for a Science Division grant. 
This work is supported in part by the National Science Foundation under grant num- 
ber PHY-9503756. 
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FIGURE CAPTIONS 



Fig. 1: The absolute square of the autocorrelation function is plotted as a function of 
time in atomic units for a particle in a two-dimensional box with fi\ = n 2 = 18, 
Li = y/3/2, and L 2 = l. 

Fig. 2: The absolute square of the autocorrelation function is plotted as a function of 
time in atomic units for a particle in a two-dimensional box with ni = n 2 = 18, 
L\ — 1, and L 2 = a/3. 
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